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Abstract. We investigate the validity of a zeroth thermodynamic law for non- 
equilibrium systems. In order to describe the thermodynamics of the glassy 
systems, it has been introduced an extra parameter, the effective temperature 
which generalizes the fluctuation-dissipation theorem (FDT) to off-equilibrium 
systems and supposedly describes thermal fluctuations around the aging state. 
In particular we analyze two coupled systems of harmonic oscillators with Monte 
Carlo dynamics. We study in detail two types of dynamics: sequential dynamics, 
where the coupling between the subsystems comes only from the Hamiltonian; 
and parallel dynamics where there is another source of coupling: the dynamics. 
We show how in the first case the effective temperatures of the two interacting 
subsystems are different asymptotically due to the smallness of the thermal 
conductivity in the aging regime. This explains why, in structural glasses, different 
interacting degrees of freedom can stay at different effective temperatures, and 
never thermalize. 



1. INTRODUCTION 

The dynamics of glassy systems has been a subject of intensive research |[. 
In the last decades there has been an increasing interest in trying to extend the 
thermodynamics ideas to systems which are far out of equilibrium. The general 
problem is to quantify non-equilibrium fluctuations within a thermodynamic approach. 
A salient feature of systems which are in equilibrium is the fact that the equilibrium 
fluctuations and the linear response functions are related by the so-called fluctuation- 
dissipation theorem (FDT) which does not hold for off-equilibrium systems. 
Several studies of spin-glass mean-field models have shown that a generalization of the 
fluctuation-dissipation theorem is possible through the definition of the "fiuctuation- 
dissipation ratio" (FDR)|, |: 

Xit,s)^^^^it>s), (1) 

ds 

which is equal to 1 in equilibrium. It turns out that the behavior of the quantity 
X{t, s) is non trivial in the limit t, s — > oo. If the lowest time s is sent to infinity the 
quantity X{t,s) becomes a non-trivial function of the autocorrelation C{t,s). This 
strong statement has been proved to hold in the framework of mean-field spin glasses 
^. Moreover, it has been recently recognized that the quantity X is generally 
related to the Parisi order parameter P{q) which appears in equilibrium studies of 
spin-glasses providing a natural link between the static and dynamical properties 

According to relation (]l|) the usual fluctuation-dissipation relation would be 
recovered if the temperature into the right hand side of (0) were T/X{t, s). This last 
ratio receives the name of effective temperature and it has been shown ||^, ^ that it 
has some properties of a macroscopic temperature. In fact, a proper thermometer 
coupled to the slow degrees of freedom can measure it. The question about the 
convenience of this temperature to describe the non-equilibrium behavior has been 
a subject of controversy in the last years |10[. While there are some evidences (not 
only theoretical but also experimental that the violation of FDT gives a 

good temperature in the thermodynamic sense, it is unclear what properties of non- 
equilibrium temperatures are common to the equilibrium ones. 

In this letter we want to analyze how effective temperatures equalize when two 
systems out of equilibrium are put in contact. In other words, we want to study if there 
exist a zeroth law for non-equilibrium systems. Let us imagine about a vitrified piece 
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of silica quenched to the room temperature. Because the glass is off-equilibrium its 
effective temperature is higher than room temperature, but when we touch the glass it 
is not hotter than the room temperature. Then we can think that some "fast" degrees 
of freedom are thermalized to the room temperature while other "slow" degrees of 
freedom remain non-thermalized. Why different interacting degrees of freedom have 
not reached thermal equilibrium for sufficient long times? We believe that the answer 
lies in the fact that the conductivity in the aging state can be extremely small, and the 
two parts of the system never thermalizep^. However, we can achieve a more deep 
understanding through a detailed analysis of an illustrative example as a previous 
stage to offer more simple and generic considerations. 

The model is a set of harmonic oscillators evolving by Monte Carlo dynamics 
introduced in jl^. The importance of this model relies on the fact that at zero 
temperature it shows typical features of glasses such as aging in correlation and 
response functions. Our interest will be in considering two coupled sets of harmonic 
oscillators. Thus, we can see how the main observables are affected by the coupling, in 
particular how the effective temperature evolves for the two sets of interacting degrees 
of freedom (i.e for the two different sets of harmonic oscillators). The interaction may 
then appear through the Hamiltonian or through the Monte Carlo dynamics itself. 
We will discover that the effective temperature for the two sets of oscillators depends 
on how the coupling is done, and we will understand why in vitreous systems different 
degrees of freedom may stay at different temperatures without thermalising at very 
long times. In this sense the utility of the extension of the zeroth thermodynamic law 
to the non-equilibrium aging state is questioned. 

The paper is organized as follows. Section II describes the main aspects as well 
as the interest of the model. Section III describes the two classes of couplings we 
have considered. Section IV discusses the results and the physical consequences of our 
work. The last section presents the conclusions. 

2. A SOLVABLE MODEL OF GLASS 

As a simple model of glass we will consider a system of uncoupled harmonic oscillators 
evolving with Monte Carlo dynamics introduced in and recently revisited in 
HQ. The Hamiltonian is: 



The low-temperature Monte Carlo dynamics of this system shows typical non- 
equilibrium features of glassy systems like aging in the correlation and response 
functions. The interest of this model is that at low temperatures the low acceptance 
rate generates entropic barriers that made the dynamics of the system extremely slow. 
The simplicity of this model makes it exactly solvable yielding a lot of results about 
the non-equilibrium behavior. 

The Monte Carlo dynamics is implemented by small movements: Xi — > Xi+ri/^/N 
where are random variables Gaussian distributed with zero average and variance . 
The move is accepted according to the Metropolis algorithm with probability W{AE) 
which satisfies detailed balance: W{AE) — W{—AE)exp{—(3AE), where AE is the 
change in the Hamiltonian. The main features of this system at zero temperature we 
are interested in are: 




(2) 
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(i) Slow decay of the energy. The evolution equation for the energy is Markovian. The 
asymptotic large-time expansion of this equation shows that the energy decays 
logarithmically E{t) ~ l/log(i). It can also be seen that the acceptance ratio 
decays faster A{t) ^ l/{t \og{t)). 

(ii) Aging in correlations and responses. The correlation function C'{t, s) is defined 
by: 

1=1 

The response function G{t, s) is calculated by applying an external field to 
the system. Then, the response function is the variation of the magnetization 
M{t) = jj- X^ili Xi{t) of the system when the field is applied: 

.>.^ ,4) 

The asymptotic scaling behavior for these two quantities is given by, 

C{t,s)^Cis,s)^ , Git,s)^Gis,s)^eit~s) , (5) 

with D{t) ~ t{\og^{t)) I and C(s,s) = ^^,G(s,s) = ^ where f{t) is a 
function that contains the whole dynamics and decays like 1/t. [|l3j. 

(iii) The effective temperature. As said in the introduction, the effective temperature 
is defined in terms of the FDR eq.dl]): 



T 



dC{t,s) 



X(t,s) G{t,s) 

As we expected, in equilibrium X{t, s) = 1 and Teg = T. Moreover, E{s) — T/2 
and the equipartition theorem is verified. In the non-equilibrium case, the effective 
temperature only depends on the lowest time s for all times s and t. This is a 
characteristic feature of this particular model that is believed to hold for structural 
glass models in the asymptotic limit s oo. 

At zero temperature when slow motion sets in, the system never reaches the 
ground state. In this aging regime the effective temperature verifies, in the limit 
s — > oo : 

TM = 2Eis) + ^ 2Eis), (7) 

/(s) ds 

Equation |^ shows how the equipartition theorem can be extended to the glassy 
regime. 

3. TWO COUPLED SYSTEMS OF OSCILLATORS 

Now we consider the case in which we couple two systems of harmonic oscillators. In 
this case it is possible to compute analytically how one system affects the other. The 
Hamiltonian is: 

2=1 i—1 i—1 

J In the original paper [ist the logarithmic corrections were estimated to be log2(t). The correct 
exponent for the logarithmic corrections was later evaluated by Th. M. Nieuwenhuizen |l5|. 
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where we take K1K2 > e^, otherwise the system has no bounded ground state. We 
define the following extensive quantities: 



where Ei and E2 are the energy of the bare systems while Q is the overlap between 
them. In this case we also consider Monte-Carlo dynamics, using the Metropolis 
algorithm for the transition probability. Although the random changes in the degrees 
of freedom Xi, yi arc defined in the same way we have explained in the previous section 
for the case of a single system, there are different ways to implement the dynamics 
depending on the updating procedure of the variables Xi,yi. In this letter we will 
analyze two different types of dynamics: 

(i) Uncoupled or sequential dynamics. In this case the two sets of variables x and 
y are sequentially updated. First the Xi variables are updated and the move 
is accepted according to the total change of energy AE = AEi — eAQ. Next, 
the variables yi are changed and the move accepted according to the energy 
change AE = AE2 — eAQ. This procedure is then iterated. In this case, the 
dynamics docs not affect simultaneously the two sets of variables but each set 
is updated independently from the other. The only coupling between the two 
sets of oscillators comes from the explicit coupling term eQ in the Hamiltonian. 
For e = the dynamics is trivial because the dynamical evolutions are that of 
two independent sets of harmonic oscillators everything reducing to the original 
model described in section II. 

(ii) Coupled or parallel dynamics. In this case the Xi,yi variables are updated in 
parallel according to the rule Xi ^ Xi + ri/y/N, yi ^ yi + Si/y/N where rj and Sj 
are random variables Gaussian distributed with zero average and variance Ai and 
A2 respectively. The transition probability for that move W{AE) is determined 
by the change in the total energy AE = AEi + AE2 — eAQ introducing, on top of 
the explicit coupling term eQ in the Hamiltonian, an additional coupling between 
the whole set of oscillators through the parallel updating dynamics. Contrarily 
to the uncoupled case, the e = case is interesting by itself because it shows how 
this kind of dynamical coupling strongly influences the glassy behavior. In fact, 
in the limiting case e = 0, there will be some changes which make the energy 
of one of the two systems increase, this change being accepted because the total 
energy will decrease. Because of that, despite of the fact that there is no direct 
coupling in the Hamiltonian the dynamics turns out to be strongly coupled. 

Our interest will focus on the behavior of two-times quantities such as correlations, 
responses and the corresponding effective temperatures. These quantities will refer to 
three classes of systems: the set of oscillators described by the x variables, the set 
of oscillators described by the y variables and the whole set of x and y variables. In 
the rest of the paper the subindex 1 will refer to quantities describing the set x of 
oscillators, the subindex 2 will refer to quantities describing the set y of oscillators 
and the subindex T will refer to quantities describing the whole set of oscillators x 
plus y. The main set of correlation and response fmictions we are interested in are: 

• Correlations. The correlation function for the sets x and y, 




(9) 



i=l 



i=l 



i=l 
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as well as the global correlation Crit, s) ~ i(Ci(i, s) + C2{t, *))• 

• Response functions. The response function for the sets x and y are defined in the 
following way. Consider two external fields hi and /12 conjugated respectively to 

W E^i ^^ and M2 = ^ 



the magnetizations Mi — jj X^ili and M2 = jf X^iLi 2/*' then: 



N ^ N 

H = -^^a;,^ + -^^yf -^(/iiXj + /i22/.0-e^Xiyj .(11) 

i—1 i—1 i i 

The responses Gi{t, s), G2{t, s) measure the change in the magnetizations Mi{t) 
and M2(t) induced by their respective conjugated fields hi and /12 applied at time 
s. These are defined by 

G,,(M) = fi^^l . (12) 



Shi^2is) 



'11.2=0 



Apart from these two response functions we may define the global response 
function GT{t,s) as the change in the global magnetization Mt — ^(Afi + M2) 
induced by a field conjugate to the total magnetization, 



GTit,s) 



• Effective temperatures. From the correlation and response functions we may define 
three effective temperatures: T^g for the system 1, T^g for system 2 and for 
the global system. These are defined as follows, 

/ aci(t,s) \ / dC2(t,s) \ / acTjt.s) \ 

We will analyze in detail the three effective temperatures for the coupled and 
the uncoupled cases, and show in which cases, depending on the dynamics, these 
effective temperatures equalize. 

4. RESULTS 

In this section we will summarize and analyze the results obtained in Ref. [p^ where 
the interested reader will find all the technical details. It can be shown that the 
equilibrium results are the expected ones. Independent of the dynamics, the effective 
temperatures are just the temperature of the bath: 

T.iff = (2i?i - eQ)/7V = T, T^^ ^ i^E2 - eQ)/N ^ T. (15) 

Where, in equilibrium, the energies of the subsystems are: 

Now we analyze the results for the effective temperatures for the two different dynamics 
defined in the off-equilibrium regime. 
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4-1- Sequential case 

The first quantity we have to focus on is the energy of each of the two subsystems. 
In this case we can find an asymptotic solution for the dynamical equations At 
first order in logarithmic corrections 1 / log(i) we find in the limit e « : 

' m{K,K2-e^)\og{t) ' ' 16{K,K2-e^)\og{t) ' 

Solving the dynamic equations for the correlations and responses one can find 
an analytic expressions for the effective temperatures. Considering the two times t, s 
both large but t — s -^i s. For a weak coupling (i.e e w 0) the value of the effective 
temperatures are, in the limit s — > oo: 

t;. . (mM)/N + o(,-) » (18) 

This yields in the s oo limit a non vanishing relative difference T^ff/T^g—l. This 
is a consequence of the fact that the two energies are different in the long-time regime 
(p^. For each subsystem the effective temperature verifies the equipartition theorem 
for long times. This means that, at any time, each subsystem can be considered as if 
it were at "quasi-equilibrium" at their corresponding effective temperature. The fact 
that the two effective temperatures are different implies that there are some degrees 
of freedom hotter than others. One can then imagine that there is always some 
kind of heat transfer or current flow going from the "hot degrees" of freedom to the 
"cold" ones. Then, one may ask why the effective temperatures do not asymptotically 
equalize. The reason is that the off-equilibrium conductivity may vanish with time 
fast enough for the heat transfer not to be able to compensate such difference ||. In 
this situation, if we now compute the total effective temperature for the whole system 
we see that in the off-equilibrium regime this temperature does not coincide with the 
sum of the energies of the systems. This is the same situation we found in equilibrium. 
If we have two systems in local equilibrium at different temperatures and put them 
in contact, the global system never verifies FDT unless the two temperatures are 
the same. In our case, we have two systems which are in "quasi-equilibrium" at two 
different effective temperatures, so the would never be the sum of the two energies 
unless the two effective temperatures T^fi,T^g were the same. In other words, two 
systems thermodynamically stable at different temperatures are not globally stable 
when put in contact. 

4-. 2. Parallel dynamics 

As in the case without coupling, the interesting dynamics is when the temperature of 
the bath is zero. In this case, the energies and the overlap decay to zero logarithmically: 

' 8{K,K2~e^)log{t) ' ' 8{K,K2-e^)log{t) ' 

Q = - ^V--iV with + (20) 

§ A detailed study can be found in [|l2|] 
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For the effective temperature up to order we may write, in the limit s — > oo (with 
^ finite): 

T^s^{2E,{s)-eQ{s))/N , T^g ^ {2E2{s) ^ eQ{s))/N, (21) 

because the asymptotic values of the Ei{s) and £'2(5) are the same the effective 
temperatures for the subsystems are also the same in the long-time limit. Note that 
the case with dynamic coupling or parallel dynamics is qualitatively different from the 
case without dynamic coupling or sequential, because now all the degrees of freedom 
are at the same effective temperature in the long-time limit. As a result, if we consider 



the global system, the total effective temperature defined in (14) is, in the limit s ~> 00 
with - finite: 

s 

Tj^ = = r,2ff = {2E'{s) ~ eQ{s))/N (22) 

where E' ~ Ei — E2 and Q are given by (pO|). This is a consequence of the fact 
that the energies of the two systems equalize due to the dynamic coupling. Now the 
dynamic coupling equalize the energies of the subsystems and the whole system has 
the same effective temperature and we can define an effective temperature for the 
global system using FDT. The situation is the same as in equilibrium systems. If we 
have two systems in equilibrium at a certain temperature T, FDT not only holds for 
each subsystem but also holds for the whole system bringing the temperature of the 
bath T. If we restrict to the case in which the coupling constant vanishes, e = 0, 
then the systems are still coupled only through the dynamics and we obtain the same 
quahtatively results: T^g = T^g = T^g = i2E'{s))/N with E'{s) w WT^^- We 
conclude that the dynamic coupling does not allow the presence of more than one 
effective temperature in the whole system because even in the absence of explicit 
coupling in the Hamiltonian, the dynamics itself makes the energies to equalize in the 
long-time limit regime. 



5. CONCLUSIONS 



In this paper we have focused our attention on the concept of the effective temperature 
defined through the FDR (|l|). The effective temperature, a parameter defined as an 
extension of FDT to the off-equilibrium regime has been introduced in the context 
of glass theory in order to understand the physics behind the dynamic behavior 
of these out-off-equilibrium systems. In this paper we hope to have clarified some 
aspects behind the physical meaning of this effective temperature. We have studied 
two types of couplings between the two subsystems of oscillators, both in an aging 
state, finding that the way we couple them is crucial for the validity of the zeroth law 
in the off-equilibrium regime to hold. The two cases we studied are the dynamically 
uncoupled or sequential case and the dynamically coupled or parallel case. In short, 
for the sequential case the coupling between the variables of the two subsystems in 
the resulting dynamics arises only through the Hamiltonian term eQ. For the parallel 
case, the variables of the two subsystems are simultaneously updated leading to further 
interaction between the two subsystems (on top of the eQ coupling term in the energy) . 

For the dynamically uncoupled or sequential case the two subsystems 
asymptotically reach different effective temperatures which never equalize. So we 
can divide the system into two parts, each part characterized by its own effective 
temperature. The reason for this behavior is that off-equilibrium thermal conductivity 
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decays very quickly to allow for an asymptotic equalization of the two effective 
temperatures. In fact, it decays as p2[ : 

Physically this means that, as time goes on, although the effective temperatures 
of the systems are different the conductivity is not high enough and they cannot 
thermalize. This seems to be the reason for the existence of different temperatures 
in real glasses, in which there are fast degrees of freedom and slow ones at different 
effective temperatures which never equalize due to the smallness of the conductivity. 
Our conclusion is that the zeroth law is probably valid but hardly effective due to the 
very small conductivity between the two subsystems in the aging state. 

We have seen that for the dynamically coupled or parallel case, the two effective 
temperatures equalize and the two subsystems are in a sort of thermal equilibrium 
between them in the aging state. As a result, all degrees of freedom have the same 
temperature. In this case, the direct coupling of the two subsystems through the 
parallel dynamics makes the conductivity much larger than in the sequential case so 
in this case a zero-th law for the aging state is effective and holds. In fact, these results 
are also valid when we consider the particular case e = in which the dynamics in 
itself is enough to equalize the effective temperatures. 

Dynamics in real structural glasses involves short scale motions of atoms and 
coupling between the different degrees of freedom occurs at the level of the energy or 
Hamiltonian and never at the level of the dynamics. Then, from these two types of 
couplings the first one is the only realistic. 

The results of this paper explain then why different degrees of freedom in 
structural glasses can stay at different effective temperatures forever. The off- 
equilibrium conductivity or heat transfer between the different degrees of freedom 
is small enough for the equalization of the effective temperatures associated to the 
different degrees to never occur. This explains why when we touch a piece of glass 
we feel it at the room temperature. To conclude, although a zero-th law for non- 
equilibrium glassy systems may hold, it is hardly effective because of the small energy 
transfer occurring between degrees of freedom at different effective temperatures. 
It would be very interesting to pursue this investigation further by studying other 
solvable examples and showing that what we have exemplified here is a generally valid 
for structural glasses as well as for other glassy systems. 
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